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Abstract 



We study the ground state properties of the double-exchange systems. The 
phase factor of the hopping amplitude arises from t2 g spin texture in two or 
more dimensions. A novel "flux" state has a lower energy than the canted 
antiferromagnetic and spiral spin states. In a certain range of hole doping, 
a phase separation occurs between the "flux" state and antiferromagnetic 
states. Constructing a trial state which provides the rigorous upper bound on 
the ground state, we show that the metallic canted antiferromagnetic state is 
not stable in the double exchange model. 
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The double-exchange model |T|-§[ has received special attention since the discovery of the 
colossal magnetoresistance (CMR) in manganites [|J. The least common model includes a 
single tight-binding band of e g electrons with the nearest neighbour hopping (t) coupled with 
the localized t 2g spins (S) by the strong Hund coupling (J#). The t 2g spins interact each other 
by the antiferromagnetic (AF) exchange interaction (J). The hopping amplitude depends 
on the relative directions of the localized magnetic moments. Anderson and Hasegawa 
derived the elements from the analysis of two sites problem. Intuitively, they are the inner 
product of the classical spin vectors on the adjacent sites. 

De Gennes obtained [^J the phase diagram as a function of temperature (T) and hole 
concentration (x). At zero temperature and x=0, the system exhibits the AF phase. The 
ferromagnetic phase away from x=0 was attributed to the double-exchange mechanism. The 
existence of the canted AF phase between them is an interesting and remarkable phenomenon 
induced by a competition between the kinetic energy of e g electrons and the AF exchange 
coupling among t 2g spins 0. The canted angle changes continuously between the two 
phases as a function of doping. These properties have been theoretically conceptualized 
as fundamental and inherent properties of the double-exchange model. Recently, Yunoki 
et al. obtained the phase diagram by an intensive numerical study. The results suggest 
the existence of a phase separation |5]]. (See also J7J.) In the other publications, it was 
claimed that one needs an additional phase factor in the hopping amplitude ||, and several 
properties associated to it have been reported |p|HTT[ . 



In this letter, we investigate the ground state properties of the double-exchange model in 
two dimension in the large limits of Jh and S. The model is reduced to a tight-binding model 
with complex hopping and with a single band. The phase factor of hopping arises from the 
t 2g spin texture in two or more dimensions || . (In one dimension, it can be always eliminated 
by a local gauge transformation and was missed in two sites problem 0.) We focus on the 
phase degrees of freedom, which makes remarkable effects, and obtain the following results: 
(i) A novel "flux" state has a lower energy than the canted AF M and spiral spin states 
|l^J13| . This state has a commensurate spin texture. We have optimized the shape of the 
density of states in addition to the band width and the exchange energy because the flux 
may induce an extension of the Brillouin zone or the Fermi surface instability (gap opening). 
The mechanism is similar to the flux phase discussed in the Hubbard model [|l4] , |i"5f and the 
generalized Peierls instability |16[ exhibited by lattice fermions in a magnetic field, (ii) In 
a certain range of x, the phase separation occurs between the "flux" and AF states. We 
propose a fundamental mechanism of the phase separation. 

Due to the above effects, it is concluded that the metallic canted AF state is unstable 
in the double-exchange model. More precisely, we explicitly construct a trial state whose 
energy is lower than that of the canted AF state. It provides the rigorous upper bound on 
the ground state energy, which is sufficient to prove the absence of the canted AF state for 
a wide parameter region. The dimensionality (two or more) plays an essential role and the 
mechanism is completely distinct from that observed in one dimension These results 

are stable under the perturbation with respect to the three dimensional layered structures 
if the inter-layer hoppings are small enough. The recent neutron scattering experiments 
18| , |T9|| and the transport properties are discussed in the light of the present results. 
The double-exchange model at finite doping is defined by 
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where q is the annihilation operator of spinless fermion, the summation is over nearest- 
neighbor sites in the square lattice, and i (or j) denotes the lattice site. The second term is 
the AF exchange coupling between t 2s localized spins, where Si is the classical spin directing 
(9i, (pi) in the polar coordinate P0| . The shape of the density of states depends on the phase 
factor in the hopping. Therefore, to study the ground state, we have to optimize the band 
width, the exchange energy, and the shape of the density of states . It should be noted 
that any analysis which assume an ordering pattern of two sublattices misses the possibility 
of the "flux" state a priori, because the effect of the phase factor originates from the global 
structure of t 2ff spins. We begin with parametrization of the phase factor in the hopping 
amplitude assuming the a/2 x a/2 structure of the t 2g spin texture. The extension to bigger 
structures is straightforward. We choose the parametrization of the t2 9 spin texture 
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where m 

(#iA,0)= 



and n 

=(0,0,0) 



are integers. The ferromagnetic and AF states are expressed by 
and (6 1 !, 6*2, 0) = (0, 7r, 0), respectively. The e g electrons acquire a flux 



whose magnitude is equivalent to the surface of the unit sphere surrounding by the four 
vectors @ in a motion around the plaquette. The Brillouin zone is extended to two sites 
per unit cell due to the phase degrees of freedom. The dispersion relation is 



E± = ±2ty A 2 (cosk x + cosk y ) 2 + B 2 (cosk x — cosk y ) 2 , 

where A=cos#icos#2 and B=sin8i sin #2- The Brillouin zone is (k x ,k y ) e + k y \ < 7r} 
n {\k x — k y \ < 7r}. The corresponding density of states is 
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where is the complete elliptic integral, a= 4t 2 (A 2 + 5 2 ), f3= 4t 2 (A 2 — B 2 ), and 7= 

64t 4 A 2 B 2 . (See Fig. |l|.) The splitting of the density of states corresponds to an extension 
of the Brilliouin zone. The formation of the (pseudo-) energy gap is crucial to stabilize the 
"flux" state against the canted AF state. The state gains the total energy optimizing the 
commensurate spin structure through the phase degrees of freedom (a Peierls phase). 

To obtain the phase diagram, we compare the total energy among the states exhibited 
by the Hamiltonian ([[]) with the parametrization (||p. The total energy is optimized in 
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the parameter space, JS 2 , 9's, and 0, setting t=l. (It should be noted that the energy is 
invariant under the global U(l) transformation with respect to 0, which plays an important 
role in the stability for the three dimensional layered structures as we will show below.) The 
canted AF state is not included in the (0) and we obtain the energy from the model ([!]) 
with <pi=0 and 9i—9 where the canted angle is specified by the 9. The results are shown 
in Fig. 0(a). (I) near half filling: For J5* 2 ~0.13, the "flux" state has a lower energy than 
the ferromagnetic, canted AF, and spiral spin states. In Fig. [5](a), we show energy of the 
states as a function of x at JS 2 =0.3. As JS 2 decreases, a transition to the ferromagnetic 
state occurs. This transition is due to a level crossing between them. In Fig. ^(b), we show 
energy of the states as a function of x at JS 2 =0.1. The configuration of spin angles is shown 
in Fig. f|. As JS 2 increases, the configuration asymptotically approaches to that of the AF 
state. (II) low doping region: A "flux" state with (#i, # 2 , 0) = (O, 0) and the canted AF 
state degenerate, but these states are unstable against the phase separation as we show 
below. 

The density of states of the ferromagnetic, canted AF, and spiral states has one van Hove 
singularity at the band center, while the "flux" state has two singularities and a pseudo-gap. 
When the Fermi level is in the pseudo-gap the "flux" state gains its energy. We have to note 
the difference of unit cell between the hopping amplitude of e g electron and the ordering 
structure of the t2 9 spins. The number of sites in their unit cells are not necessarily the 
same. For example, the unit cell of the ti g spins has two sites for the canted AF state, and 
a lot of sites for the spiral spin states. On the other hand, the unit cell of e g electrons has 
only one site for the canted AF, spiral, ferromagnetic spin states. 

So far, we restricted the spin structure to be a/2 x a/2. The optimized state has a staggered 
"flux" structure because the total phase around a plaquette by (2) is opposite in sign for 
two neighbor plaquettes. The state is, so called, the staggered flux state and is distinct 
from the flux state in the usual sense. Therefore, we use the double quotation to distinguish 
it. Generalizing the mechanism of the a/2xa/2 "flux" state, a flux state characterized by a 
bigger unit cell structure could have a lower energy than the canted AF state at an arbitrary 
rational filling. Because, if we assume a unit cell structure axb=q (a, b>y2), the structure 
allows to suppress the density of states at x=p/q where p and q are coprimers. Especially, 
when we restrict the discussion inside the "Flux" region in Fig. 0(a), the ax b flux state could 
gain some energy with respect to the \/2x\/2 "flux" state except at half-filling. Although 
we have not explored in such bigger structures in detail, it would be interesting to study 
the relation among the flux (or "flux") state, the rational hole concentration, and a bigger 
structure of the unit cell along the line of the present paper. 

Assuming the ordering of the t 2g spins, the phase diagram is obtained (Fig. 0(a)). How- 
ever, the canted AF phase is unstable under phase separation. The phase diagram is shown 
in Fig. 0(b). As shown in Figs. 0, the system can gain the energy by mixing the AF and 



the state whose energy level touches the tangent line from the AF point, i.e. x=0 |22 
For JS 2 ~ 0.23, it exhibits the phase separation with the AF and "flux" spin states. For 
JS 2 ~ 0.17, the system exhibits the phase separation between the AF and ferromagnetic 
states. A fraction of the pure canted AF state survives in Fig. 0(b) (iii), but the state de- 
generates to a "flux" state with (9%, 6*2, 0) = (O, 9, 0) there. The nature of the boundaries 
between the phase separation and the ordered state is expected to be of the first order. The 
existence of the "flux" state is crucial to the absence of the canted AF state. Otherwise the 
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canted AF state still remains instead of the region of the "flux" state in Fig. 0(a), and the 
phase separation between AF and canted AF occurs in the region JS 2 ~ 0.23. Although 
the y/2xy/2 "flux" state is a trial state in the sense of the variational argument, it provides 
the rigorous upper bound on the ground state energy. Therefore, this is sufficient to prove 
the absence of the canted AF state in the wide parameter region. (It is possible that the 
true ground state is a "flux" state with bigger structures or another state and the phase 
separation between the state and AF state occurs.) 

The "flux" state is stable against the perturbations with respect to three dimensional 
layered structures. This is due to the following reasons: (i) the structure of the van Hove 
singularity in each layer is stable if the inter-layer hoppings are small enough, (ii) Each layer 
has the global U(l) symmetry with respect to 0, and can adjust the AF exchange energy 
taking different 0's. The quantitative study is a problem in future. 

Experimentally, the "flux" state and the phase separation between the state and AF 
state are particularly interesting. The spin configuration of the "flux" state is unusual |13| . 
In the neutron scattering experiments, this will be detected as three spots, (vr/2, vr/2), (0, it), 
and (n, 0) corresponding a complex folding of the Brillouin zone. A remnant of the phase 
separation between the "flux" and AF states may be detected by a crossover of the intensities 
as a function of x, because the intensities of each component can be evaluated independently. 
However, one cannot distinguish the phase separation from the other phenomena only from 



the neutron scattering (see |T8|,pJ|,|23| ) . The ax b flux state could gain some energy with 



respect to the \/2xy/2 "flux" state in the "Flux" region in Fig. |](a) or to the canted AF 
in the other region. It is expected that the stabilization occurs with an energy gap due 
to a similar mechanism to ]TJ| or due to the mechanism [17|. These conditions meet the 
quantization of the Hall conductance P^LKH . 

In summary, we have investigated the universal properties of the double-exchange mech- 
anism at zero temperature. The "flux" state has a lower energy than the canted AF and 
spiral states. The phase degrees of freedom and the dimensionality (two or more) play an 
essential role. It is concluded that the metallic canted AF state is unstable in the double 
exchange model within the analysis of the least common model, i.e. without any assistance 
of other effects like polaron [0], orbitals, etc. This is one of few examples of the "flux" state 
26| realized in non-singular model. 
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FIGURES 



FIG. 1. Densities of states in two dimension. The solid line is for the "flux" state. The broken 
line is for the ferromagnetic, canted AF, and spiral spin states. 

FIG. 2. (a) Phase diagram of the double exchange model, which exhibits the "flux" ("Flux") 
and ferromagnetic (F). The canted AF state and "flux" state with (0i,02,4>)=(O,9,4>) degenerate 
(CAF). (b) Phase diagram after taking into account the phase separation, which exhibits the 
phase separation between AF and "flux" (AF+'Flux"), AF and ferromagnetic (AF+F), AF and 
the canted AF (AF+CAF), "flux" and Ferromagnetic (i), "flux" and canted AF (ii) spin states. In 
(hi) the canted AF and "flux" states degenerate. 

FIG. 3. Energy of the "flux" state (i), canted AF state (ii), ferromagnetic states (hi) as a 
function of x at JS* 2 =0.3 (a) and JS 2 =0.1 (b). The tangent line (iv) from x=0 is an energy of 
phase separation. 



FIG. 4. Energy dependence of O's as a function of J at x=0.4. 
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